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By using Density Matrix Renormalization Group (DMRG) technique 
we study the ID extended anisotropic Heisenberg model. We find that 
starting from the ferromagnetic phase, the system undergoes two quantum 
phase transitions (QPTs) induced by frustration. By increasing the next- 
nearest-neighbor (NNN) interaction, the ground state of the system changes 
smoothly from a completely polarized state to a NNN correlated one. On 
the contrary, letting the in-plane interaction to be greater than the out-of- 
plane one, the ground state changes abruptly. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.30.Kz, 75.10.Pq, 75.40.Mg 

1. Introduction 

QPTs, contrarily to ordinary ones (at T ^ 0), are characterized by a 
change of ground state on continuous varying the external conditions [1]. 
Such phase transitions are often caused by frustration. One-dimensional 
quantum spin systems offer a wide variety of examples of QPT. Among 
them, the Heisenberg-like models are the most important. While ID Heisen- 
berg model with only nearest-neighbor interaction possesses a well studied 
analytic solution by Bethe anzatz [2], the addition of NNN term suffices to 
break the integrability. In such extended system a whole set of new phenom- 
ena occurs. The physics of this model has been already largely studied by 
means of many field-theoretical as well as numerical methods [3j , primarily 
in the proximity of the antiferromagnetic phase. 

However, the ferromagnet is also one of the ground states of the model 
and it is primarily promoted by the diagonal part of the Hamiltonian. If 
other terms in the Hamiltonian induce frustration, then eventually the fer- 
romagnet will "melt", and the system will undergo a QPT towards some 
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new phase. In the present work we investigate such melting driven by two 
distinct frustrating factors: i) the anisotropy of the nearest neighbors inter- 
action term and ii) the antiferromagnetic NNN interaction. 

2. Model and Method 

We study the one-dimensional anisotropic extended Heisenberg model 
with NNN interaction: 

H = -J.Y. ^tsu, + J± Y^^SfSf^, + sfsl,) + J'Y, S^Si.+2. (1) 

i i i 

We find the ground state of the Hamiltonian ([T|) numerically by means of the 
DMRG \4l technique on a chain with 100 sites, subject to open boundary 
conditions. We use Jz > 0, which corresponds to ferromagnetic coupling. 
In the present work we calculate the Fourier transform of the correlation 
functions, both in-plane {{S^{k)S^{k + q))) and out-of-plane {{S^{k)S^{k + 
q)))- Due to the large system size, the above correlation functions do not 
depend on k, but only on the transferred momentum q. The starting point 
for our analysis is the doubly degenerate completely polarized state with all 
spins either "up" or "down", located on the phase diagram of the system 
within the limits < J± < Jz and J' < O.SlJz, as found in Ref. [3]. 

3. Frustration due to J' 

Up to the values of J' 0.31J^ the system remains totally polarized. 
In the momentum space this corresponds to a 5— function peak at g = 
in z— channel and a constant value of 1/4 in the perpendicular one. On 
the other hand, in the limit J' — > oo the system behaves as a couple of 
non-interacting antiferromagnetic Heisenberg chains with nearest-neighbor 
interaction. It is clear that in the case of J' ^ oo the system will exhibit 
short-range AF correlations with exponential decay in real space, which 
corresponds, in momentum space to a couple of Lorentzian peaks centered 
at q = ±Tr/2. 

Smooth change of J' causes a smooth rearrangement of the ground state. 
Namely, for J' > 0.317^ the system abandons the totally polarized state, 
although conserving much of the magnetization along the z— axis (see the 
peaks at g = on Fig. [1^) for J' > O.SIJ^). At the same time a short- 
ranged incommensurate phase appears: two strong Lorentzian-type peaks 
develop in the vicinity of q = ±7r/2, the larger is J' the closer are the peaks 
to ±tt/2. 

On the contrary, in the perpendicular channel we have perfect commen- 
surability (at least for the values of J± < 0.5 Jz) as shown on Fig. Wp)- As 
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in the case of the z-axis, the in-plane correlations appear to decay exponen- 
tially with the correlation length depending on both J_\_ and J' . 

4. Frustration due to Jj_ 

Contrarily to the smooth change of the properties in the previous case, 
frustration caused by J_l changes the phase abruptly. The spins, totally 
polarized along z direction for J±_ < J^, become antiferromagnetically or- 
dered in the x — y plane (see Fig. Hh))- The remaining component of 
the spins remain uncorrelated already at few lattice spacings. In the limit 
J_L — > oo we expect the system to become the XX model, with an in-plane 
AF ground state. As a matter of fact, at J± = 1.2 most of the correlations 
are concentrated in the x — y plane. The scaling behaviour in the x — y plane 
is once again exponential (Lorentzian in fc-space, see Fig. [iji))- The value 
of J' does not qualitatively change this picture within the range considered 
(0 < J' < O.3J2). Beyond this range the two phases described above, should 
merge at some transition line, situated somewhere in the region J' > O.SIJ^, 
0.5 Jz < J±0.9 < Jz. 

5. Summary 

In conclusion, by using a high-precision numeric tool (DMRG) we have 
demonstrated how different types of frustration can lead to QPTs in the ID 
extended anisotropic Heisenberg model. By considering a chain with 100 
sites we have been able to observe the exponential decay of the correlations 
and obtain a high resolution picture of the static spin-spin correlations in 
momentum space. In contrast to the frustration induced by J_l, the one 
due to J' leads to a rather smooth rearrangement of the ground state. A 
more detailed study of the two phases described in the present work and 
the transition between them is in progress and will be reported elsewhere. 
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